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ABSTRACT 

A rocket thrust chamber assembly was analyzed to determine 
the exit conditions that result when heat is withdrawn from the combus- 
tion chamber and simultaneously added to. the gas flow in the nozzle. 

The analysis was performed for three cases. The first case considers 
the flow of an inert gas in the nozzle , the second considers the equi- . 
librium flow of a dissociating diatomic gas, and the third case considers 
the addition of heat to the dissociating diatomic gas flow in the nozzle 
without the removal of heat from the combustion chamber. 

Results of the inert gas analysis show that transferring heat 
from the combustion chamber to the nozzle causes a decrease in both exit 
velocity and specific impulse but results in an increase in exit pres- 

; '•« '•> • ' ■ J ..." I ' Jr.' 

sure and exit temperature. The dissociating diatomic gas behaves in a 
manner similar to the inert gas, but the magnitude of the decrease is 

01 f , w . '* > i. 

not as great. This damping effect is caused by dissociation of the 
diatomic gas. 

When heat is applied to the nozzle., but not removed from the' 
chamber, the dissociating diatomic gas shows an increase in exit vel- 
ocity, exit pressure, exit temperature, and specific impulse. 
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INTRODUCTION 


Many authors have investigated various aspects of gas flows in 
nozzles, ranging from determining the effects on performance of employ- 
ing different nozzle designs to using various combinations of fuels and 
oxidizers. Recent analytical investigations conducted by Michael Smith 
(ref. 1) have been concerned with external heat addition to the nozzle 
exhaust gases and sudden freezing of chemical and vibrational rate 
processes during flow in de Laval nozzles. 

The author of this paper has extended Smith’s work to the case 

where the quantity of heat added to the nozzle is simultaneously with- 

drawn from the combustion chamber. The analysis is extended to cover 
two cases: the first case considers the flow of an inert gas in the 

nozzle, and the second .case considers the equilibrium flow of a dis- 
sociating diatomic gas. In both cases it is assumed that the gas flows 
isentropically from the combustion chamber to a position immediately 
ahead of the heat addition plane and also from immediately behind the 

plane to the nozzle exit. If the zone of heat addition is infinitesi- 

mally thin, then constant-area heat addition equations can be used to 
represent the flow at the heat addition plane. The combustion chamber 
is treated as a large reservoir of stationary gases with heat being 
removed. The appropriate equations for each flow regime are derived 
and a check calculation of the computer solution is made. In addition, 
the case of external heat addition to a dissociating diatomic gas 
without heat removal from the combustion chamber is investigated. 
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One possible application of this analysis can be found in 
nuclear rocket engines where a large amount of thermal energy is avail- 
able. However, since the combustion chamber and nozzle material have 
certain temperature and pressure limitations, only a limited amount of 
energy can be added to the combustion chamber. Excess energy can be 
added to the flow at a plane in the nozzle at which the temperature and 
pressure are lower. The uniform addition of heat at a particular cross' 
section of the nozzle presents a major difficulty if the flow is not to 
be disturbed mechanically. One means of accomplishing this energy 
transfer is by thermal radiation; however, this thesis does not inves- 
tigate methods of heat, transfer. Uniform heat addition in an infini- 
tesimally thin plane is simply assumed and an analysis is made based on 


this assumption. 



CHAPTER I 


GENERAL DISCUSSION 

The rocket engine thrust assembly used in this paper consists 
of a combustion chamber and a nozzle. The combustion chamber is consid- 
ered as that portion of the overall assembly where the majority of the 
conversion from chemical to thermal energy occurs and where the gas 
speed is much less than its speed in the nozzle. The nozzle is, in turn, 
considered as that portion of the assembly where the majority of the 
conversion from thermal to kinetic energy occurs. A profile view of 
this configuration for a liquid rocket engine thrust assembly is shown 
in figure 1. 

Penner (ref. 2), derived the relations for the thrust and 
specific impulse of a rocket motor in the following manner: from the 

equation for conservation of momentum, the thrust on^ a rocket motor is 
given by 


F = 


(l + cos a) 


Hr v 

ex 

2 g 


+ 




(i) 


where 

a - half of the divergence angle of the nozzle 
W = the rate of total weight flow 


1 



2 


V = linear flow velocity at the nozzle exit position 

GX 

g = acceleration due to gravity 

P ex = absolute pressure at the nozzle exit 

= ambient atmospheric pressure 

A© x - cross-sectional area at the nozzle exit position 

If a is sufficiently small, equation (l) can be approximated by 


F = 


W V 


ex 


g 




( 2 ) 


The most generally accepted measure of the performance of a 

rocket motor is its specific impulse I , which is defined as the 

sp 

ratio of thrust to propellant weight flow rate. Thus, when equation (2) 
is divided by W and a is small, the specific impulse equation for a 
rocket motor is found to be 


I 

sp 


F_\x + 

W g 



( 3 ) 


Maximum thrust, and thus maximum specific impulse for a 

constant W, is obtained for ideal conditions in which the exit pressure 

P equals the external pressure P • If P = P , 
ex * am ex am 


then 



3 


equation (3) reduces to 


I 

sp 



(10 


A nozzle expansion for which P = P is sometimes referred to as 

ex am 

"being ideally expanded. 

In order to utilize equations (3) and (4) for practical eval- 
uation of motor performance, an ideal motor can be hypotesized; such a 
motor satisfies the following assumptions: 

1. Thermodynamic equilibrium is reached in the combustion 
chamber after the original adiabatic reaction. 

2. Expansion of the combustion products through the de Laval 
nozzle is adiabatic. 

3. The products of combustion behave as ideal gases. 

4. The expansion process is considered to involve one- 
dimensional flow (parallel to the nozzle axis) of non- 
viscous ideal gases. 

5. The velocity of the gases at the nozzle entrance position 
is negligibly small compared to the velocity at the exit 
position. 

This Idealized rocket motor is used as a standard of comparison for two 
specialized thrust chambers investigated in this thesis. For a more 
detailed development of the subject matter presented above, the reader 
should consult Penner (ref. 2), Barrere (ref. 3), Hesse (ref. 4) , or 
other authors of texts concerning rocket propulsion systems. 

The purpose of this thesis is to examine, for two specific 
cases, the theoretical effects on exit velocity and specific impulse 



due to a transfer of thermal energy. For this first case (chapter II), 
the effects on exit velocity and specific impulse will he determined 
when heat is taken out of the combustion chamber and the same amount 
of heat is added simultaneously to the flowing gas behind the nozzle 
throat for a steady state inert gas system. Previous work in this area 
by Michael L. Smith (ref. l) involved heat addition to the gas in the 
nozzle. Results from this work indicated that the addition of heat to 
the flow of an inert gas in a de Laval nozzle always results in larger 
values for specific impulse and exit pressure. On the other hand, the 
kinetic energy of the exhaust gases may or may not be increased, depend- 
ing on where heat is added to the flow. Heat addition ahead of a 
particular area ratio, which Smith found to be approximately half the 
exit area ratio, resulted in an increase in exit velocity. Smith also 
found that heat added downstream of the 11 critical" area ratio decreased 
the exit velocity. Thus, there exists a critical area ratio at which 
heat addition yields an exit velocity identical to the exit velocity 
with no heat addition. As the location of the heat addition plane moves 
downstream from near the throat to the exit area, Smith’s results showed 
that the specific impulse ratio and exit velocity ratio both decreased, 
whereas the exit pressure ratio increased. Chapter II is thus an exten- 
sion of Smith’s work to the case where the external heat added to the 
gas flow comes from the reservoir (or combustion chamber) of the rocket 


motor. 



In the second case (chapter III), the gas employed is 
hydrogen in equilibrium flow instead of an inert gas. The gas is 
assumed to be in chemical, vibrational, and rotational equilibrium 
throughout the nozzle. Other assumptions are the following: 

1. Ionization is considered to be negligible. 

2. The flow is choked (i.e., the Mach number is unity at the 
nozzle throat). 

3* No flow separation occurs in the system. 

The calculated exit velocity and specific impulse values of the inert 
gas case and the calculated exit velocity of the real gas case are com- ' 
pared with their counterparts for similar initial conditions for an 
idealized thrust chamber. ■ The idealized system is one in which the 
assumptions listed for a hypothesized ideal motor are valid and the 
flow is ideally expanded (i.e., exit pressure equals ambient pressure). 
The ambient pressure for both cases is constant for all calculations. 



CHAPTER II 


INERT GAS CASE 

A. General Considerations 

This chapter considers the equilibrium flow of an inert , ideal 
gas in a thrust chamber. If the chamber pressure is greater than the 
critical chamber pressure needed to achieve sonic conditions at the 
throat , the flow upstream of the throat is subsonic and downstream of 
the throat is supersonic. The specific impulse of this system is given 

by equation (it) for an ideal expansion process (i.e., P = P ). 

6x am 

A predetermined amount of thermal energy q per unit mass 
flowing through the nozzle is removed from the reservoir and simultane- 
ously added to the flow in the nozzle at a specified location. Removing 
heat from the reservoir results in a decrease in the temperature of the 
gas. The chamber pressure, meanwhile, is assumed to remain constant. 

The next assumption is that the gas flows isentropically from the res- 
ervoir through the throat and up to the plane where the heat is being 
added. At the heat addition plane, the heat .is assumed to be added to 
the flow in an infinitesimally thin region and, as a result, the thermo- 
dynamic state variables and the gas velocity are discontinuous across 
the plane. Immediately behind and downstream of the heat-addition plane, 
complete equilibrium is assumed and the gas expands isentropically to 
the nozzle exit. 


6 
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The objective is to determine the conditions at the nozzle 

exit (i.e., V , T , P ,M ,1 ) when the initial reservoir condi- 

ex ex ex ex sp 

tions (P^, T^), the heat exchange parameter q, the ratio of the nozzle 

exit area to the throat area (A /A ) , and the area ratio of the heat 

ex 

** 

addition plane (A^/A ) are known. The results will he compared with 

those for an identical system hut with zero heat exchange (i.e., one 
with an isentropic, adiabatic, and ideal expansion process). Note that 
for a fixed nozzle configuration, the exit pressure P will not neces- 

6X 

sarily he equal to the ambient pressure P flm for the case of heat addi- 
tion. Thus, when q is non-zero, the effect of the pressure difference 
term in equation (3) must he included when calculating the. specific 
impulse. 

B. Mathematical Development 

As seen in figure 2, the flow in the nozzle can he divided 
into three discrete regions: isentropic flow from the entrance of the 

nozzle to the area ratio where heat is being. added to the flow, flow 
with heat addition across the discontinuity in the nozzle, and isentropic 
flow .from immediately downstream of the discontinuity to the nozzle exit. 
The necessary equations for each of these regions will he derived, as 
well as the appropriate equations for heat removed from a large reser- 
voir of stationary gas at constant pressure. 
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1. Basic Equations for the Reservoir 

The equation for conservation of energy equation for the gas 
in the reservoir may be written as 

h i " h f = q » (5) 

or, using the caloric equation of state for perfect gas, as 

Cp (t ± - T f ) = i ,■ (6) 

where is the temperature of the gas in the reservoir after the 

/V 

thermal energy q has been removed and is the initial reservoir 

temperature. Since the heat exchange parameter q is defined as 


1 = c"V ’ (7) 

p i 


equation (6) may be rearranged to yield 


T f = T. (1 - q) (8) 


For the type of reservoir specified, and represent stagnation 

temperatures. As mentioned earlier, the reservoir pressure will be 


assumed to remain constant (i.e., P. = P_). 

’if 
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2. Basic equations for Isentropic Flow Regimes of the Nozzle 
The assumption of isentropic flow and the requirement that 
mass, momentum, and energy he conserved can he used to derive the follow- 
ing relations for the isentropic flow regimes in the nozzle, see Hesse 
(ref. k): 


P_ 

* 

P 



Y f 1 
2 





(9) 


Y + I 

T_ 2 

* y - 1 2 

T ! + M 


( 10 ) 


A_ 

* 

A 


1 

M 



j 


( 11 ) 


and 



where M, the Mach number, is defined hy 


(12) 


M = V/a 


(13) 
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and the speed of sound a for a perfect gas is 

a = (g y RT) 1/2 • (14) 

In this nozzle configuration, there exist two isentropic flow 

regimes. For the flow between the chamber and the heat addition plane 

* 

located downstream of the throat, A is the throat area. For the flow 

* 

behind the heat addition plane, A is an imaginary throat area at 
which the Mach number would be unity for the new flow. 

3. Basic Equations for Rayleigh Flow in Constant Area Ducts 
The assumption that the flow properties are discontinuous 
across the heat addition plane is equivalent to assuming that this re- 
gion is infinitesimally thin. This observation leads to the conclusion 
that the equations developed for the flow of an inviscid, inert gas in 
a constant area duct with heat addition, can be used for this region. 
This type of flow is often referred to as a Rayleigh process. The 
Rayleigh equations used to relate flow parameters across the plane of 
heat addition are, see Hesse (ref. U): 


P _ y + 1 

P 1 + Y M 2 


T 

** 

T 


M 


Y + 1 
1 + Y M 2 


(15) 


(16) 
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and 



T 

T 

where' — is the ratio of total temperatures. The above 
T 

T 

were derived from the continuity equation 


pVA = Constant 


the equation of state for a perfect gas 


P = pRT 


the conservation of momentum equation 


2 ' 2 
p + D V = P + 0 V 
1 P 1 1 2 p 2 y 2 


and the conservation of energy equation 


V 2 V 2 

Vi * -i- - c p t 2 + + 


Equation (2l) may be rewritten as 
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(IT) 


relations 


(18) 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 
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where the total temperature is defined by 

2 

t t = T + rr <23) 

P 

Equations (5) through (23) can he used to solve for values of. all 
the flow parameters at the nozzle exit. 

C.- Example Check Calculations and Comparison With Computer Results 
Several check calculations were performed for the computer 
program with the result that the answers calculated by hand compared 
very favorably with those calculated by the computer. One of these 
check calculations was performed for the following set of conditions: 

A 

-f 2 - = 20.0 

A 

^1 

= 10.0 

A 

T. = k000° R 

l 

= 30 atm , 

y = 1.U0 
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and 


q. = 0.25 


The following is the method used to perform the check calculation for 
the above conditions . 

A series of calculations is first performed that yields values 
for the reference conditions (i.e., for the case of q = 0). 


ft ' ! 5 •<: y + Is, 

L.iiiM 2 ' 2(Y - 1> 

1 .2 ex 0 


ex _ 1 

# = M 

A ex 


Y + 1 

2 


= 20 


The preceding equation -is satisfied when; M ■ 4.7245. The values of 


ex 


T and P are calculated by 
ex ex 

o o 


T = T. 
ex l 
o 


1 + 1 ~ ^ M 2 
2 ex 


i= 4000 (0.18303) = 732.131° R, 


and 


P = P. 
ex i 
o 


+ 1 ~ ~ L M 2 
2 ex 


Y - 1 


30 


381.1966 


= 0.0787 atm 



lb 


Equation (22) is now used in the following manner for the case 
of heat addition. 


T f = T. (1 - q) = 1+000 (0.75) = 3000° R 
P^ = P^ = 30 atm. 


For an area ratio of heat addition equal to 10, the following equation 
is solved for the value of 



1 + 




Y + 1 
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This equation gives a value for of 3-92313. The values of T^ 

and P^ are calculated by 


T = T 
1 f 


1 + 


- 1 ? 
-r — M„ 


= 3000 (0.21+52) = 735-621° R 


and 


P = P 
1 f 


y - 1 


30 


1 + x t Am i 


136-9735 


= 0.21902 atm 
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Equation (l8) gives the following two relations 



Since T = T. and 1L, = T„, equation (22) also implies that 

±2 i ^ 


T = T (1 - q) 
1 2 - 


or 




+ q T. 

^ l 





Thus 



l6 


Mg is then found by solving the following equation 



For the values of M^, q, It , and T T previously determined, this 
equation yields Mg = 2.0163. The values of Tg and Pg are calcu- 
lated by 


T = T 

2 V 


(i *^ 4 ) 


= 2206.166 °R 


and 


P 


2 




= O.T3T9T6 atm 


The pressure relation was found by dividing 
equation (15). 



as given by 




* 


A 



10 

20 



1 + 


1 + 



1 


1 



V + 1 

2(y - 1) 
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Solving for M yields M = 2.77573. The values of T 

GX GX GX 

P are found hy 
ex 


and 


T = T c 
ex 2 


1 + 1 T iM L 


= 1574.22° R 


> 



1 + 


1 + 




ex 


X _ 

y - 1 

= 0.22648 atm 


Forming the ratios of interest gives the following: 


T 

ex 

o 


1574.22 ■ 
732.131 


2.150189 


M 

ex = 2.77573 

M 4.72415 
ex 

o 


0.587562 


P 

ex 


P 

ex 

o 


. 2261+8 

.078699 


2.87777 


and 
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and 


V 

ex 


V 

ex 

o 


M 


ex 


M 


ex 


y g RT 
Y g RT 


ex 


ex 


o 


1/2 


0 . 86157 


9 


I V /p 

. gP- = ex_ ex 

I V l P 

sp ex^ \ ex 

r o o \ o 


1 



0.92167 . 


For the above initial conditions, the computer solution gave the follow- 
ing values for the ratios of interest: 


T 

~~ = 2.1512 

ex 

o 


M 

iT^= 0.587^0 

ex 
- o 


P 

p ^- = 2.88323 

ex 

o 


V 

y 22 - = 0.86153 

ex 

o 
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and 



0.92177 . 


D. Results and Conclusions 

The results of applying the mathematical relations developed 
in section B and used in the formulation of' the computer solution are 
presented in figures 4 through 13. 

Figure " k shows the effect of the location of the 'heat addition 
plane on the exit velocity ratio for several different values of the 
heat exchange parameter q. Figure 5 is a plot of heat addition loca- 
tion versus exit velocity ratio for three values of the exit area ratio 
A /A . Figure 6 shows the effect -of the heat addition plane location 

GX 

on the exit pressure ratio for several values Of q while figure 7 
shows the effect of different exit area ratios. Figures 8 and 9 are 
similar to figures 4 and 5 except that the specific impulse ratio is 
plotted instead of the exit velocity ratio. Similarly, figures 10 and. 
11 show the exit temperature ratio and figures 12 and 13 show the exit 
Mach number ratio. 

Figures 4, 8, and 12 show that the exit velocity ratio, speci- 
fic impulse ratio, and exit Mach number ratio decrease as the heat ex- 
change parameter q increases and also as the heat addition plane area 
ratio increases. The three ratios approach unity for low values of q 
and when the area ratio of heat addition is near unity. Figures 5> 
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9 9 and 13 show that this trend is found for all three exit area ratios 
considered. As for the exit pressure ratio and exit temperature ratio, 
figures 6 and 10 show that these ratios increase as either q or the 
area ratio of heat addition increases. The ratios are again closest to 
unity for the heat addition plane located near the throat of the nozzle. 
Figures 7 and 11 show that for the three exit area ratios considered, 
the curves are homothetic. That is, the size of the nozzle has very 
little effect on the exit parameters. Changing the value of the heat 
exchange parameter causes the greatest variation in the exit parameters 
of interest . 

Note that there. is no "crossover point" as was found hy 
Smith in his work (ref. l) on heat addition to the nozzle (i.e., there 
is no location at which heat addition yields an exit velocity identical 
to the exit velocity for the case of no heat addition) . To transfer 
heat from the combustion chamber or reservoir and add the thermal energy 
somewhere downstream of the throat, it is best to transfer as little 
thermal energy as possible and then introduce the heat back into the 
gas flow as close to the throat as practical, providing the desired 
exit parameter to be maximized is the exit velocity or the specific 
impulse. If, instead, the exit pressure or exit temperature is to be 
maximized, then the amount of heat transferred should be large and 
should be added to the gas at the nozzle exit. 



CHAPTER III 


REACTING GAS CASE 

A. General Considerations 

This chapter considers the equilibrium flow of hydrogen gas in 
a rocket engine thrust chamber assembly.- If the gas is introduced into 
the combustion chamber at a sufficiently high temperature T^, dissocia- 
tion of the molecular species will occur in accordance with the 

following endothermic reaction: 

L + H 2 * 2H + L , (24) 

where L is a third body (either H or Hg). The dissociation reaction 
will be in chemical equilibrium with the reverse exothermic recombina- 
tion reaction 

L + 2H — Hg + L , (25) 

if sufficiently low flow velocities are maintained in tile combustion 
chamber! For equilibrium, the relative concentrations of -the atomic and 
diatomic species are dependent on the pressure P . and temperature 
existing in the chamber. Ionization of the hydrogen atom, which becomes 
appreciable for temperature greater than 10 000° K, is assumed to be 
negligible for the pressure-temperature combinations considered in this 
thesis. 


21 



22 


The velocity of the dissociated gas will increase as the gas 
enters the convergent nozzle section and begins its expansion to the 
nozzle exit. The gas will be considered to be in equilibrium flow from 
the combustion chamber to immediately in front of the heat addition 
plane {where a known amount of thermal energy is added) and also from 
immediately behind the plane of heat addition to the nozzle exit. The 
terminology "equilibrium flow" refers to the condition that the chemical 
reactions as represented by equations ( 2 *+) and (25), in addition to 
other relaxation rate processes (e.g., vibrational, rotational) proceed 
at a sufficiently rapid rate to adjust to the equilibrium conditions 
corresponding to the local pressure P and temperature T. in the 
nozzle . 

The flow system of the thrust chamber assembly consists of 
hydrogen gas entering the combustion chamber at pressure P_^ and tem- 
perature T.. A predetermined amount of thermal energy is then with- 
drawn from the combustion chamber and simultaneously added to the 
flowing gas at a known plane downstream of the throat. It is assumed 
that the gas in the chamber immediately attains a new equilibrium 
temperature T^ and pressure after the heat removal. The pressure 

P^ is further assumed to be equal to the initial pressure P^. The 
dissociated gas then flows from the combustion chamber through the 
nozzle throat to the plane of heat addition. Throughout this portion 
of the flow, the molecular weight W, the specific heat at constant 
pressure C^, and the ratio of specific heats y of the equilibrium 
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gas mixture vary during the expansion due to the shift in equilibrium 
concentrations of the species H and as well as to the shift in 

internal energy distribution of the gas. At the heat addition plane, 
the thermal energy extracted from the gas mixture in the reservoir is 
added to the gas flowing in the nozzle. It is again assumed that the 
gas flow reaches a new equilibrium temperature T^ and pressure 
immediately after heat addition. The resulting gas mixture then expands 
isentropically to the nozzle exit. 

The major difference between the reacting-gas case and the 
inert-gas case is that, in the former, there is an added internal 
"heat addition" effect caused by recombination of the hydrogen atoms. 

B. Mathematical Development 

As seen in figure li, the flow in the nozzle can be divided 
into three different flow regimes. These consist of isentropic flow 
from the- nozzle entrance to a position immediately in front of the 
heat addition plane, flow with heat addition .across the discontinuity 
in the nozzle, and isentropic flow from immediately downstream of the 
plane to the exit. The problem of heat removal from a large reservoir 
of gases at constant pressure and zero velocity (i.e., the combustion 
chamber) is also considered. 

1. Calculations of Thermodynamic Properties at a Given Point 
in the Flow •' * ' 

If the temperature and pressure at a point in the equilibrium 

gas flow are specified, it is possible to find K^, Cp , Cp , 

H H2 
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Sjj, and "by using tables such as the JMAF Thermochemical Tables 

(ref. 5) or those furnished by the National Bureau of Standards. Once 
the values of the thermodynamic properties are known for molecular and 
atomic hydrogen at the specified conditions, the mole fraction (i.e., 
and and the thermodynamic properties of the gas mixture can be 

calculated by using the following sequence of operations as described 
by Penner (ref. 2). The degree of dissociation B may be calculated 
using the equation 


B = 



(26) 


Once g has been calculated, the mole fractions are given by 

^ = m ’ 

X H2 = i + e 9 ^ 28 .^ 


or 


^2 = 1 - X H 


( 29 ) 


Now that the mole fractions have been calculated, it is possible to 
calculate the entropy, enthalpy, and coefficient of specific heat at 
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constant pressure for the gas mixture using relations derived by Penner 
(ref. 2) for a mixture of ideal gases: 


and 


N 


w = £ X 3 W J 


j=l ’ 


(30) 


K 


N 


S = 


£ *•. ;i .i - 3 ( 1:1 p * £ . y ta 'v 


j=i 


j=i 


(3i) 


H " .2. X j ( AH f + H j " H j, 298 .i6°k) » 
j=i 7 


(32) 


N 

C = V X. C 
p o p 

j=l J 


(33) 


The ratio of specific heats for the specific conditions can then tye cal- 
culated by the following equation 


Y = 



(3b) 


2. Equations for Heat Removal From the Reservoir 
Assuming that the reservoir is initially at equilibrium with 
a temperature and P^ , the enthalpy of the gas mixture and any of 

the other desired thermodynamic properties can be calculated by using 
the above procedure (Section B, part l). Next is the problem of 
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determining the new equilibrium reservoir conditions after a given 
amount of thermal energy Q, has been removed from the reservoir. 

Starting with the conservation of energy equation and making two assump- 
tions that the pressure of -the reservoir remains constant and that the 
velocity of the gas mixture in the reservoir is zero, the equilibrium 
temperature in the reservoir after heat removal, T^, is found by a 
trail and error solution. The conservation of energy equation with the 
above assumptions can be written as 

dh = d 4 , (35) 

which, when expanded mathematically, becomes 

h ± - h f = 4 (36) 

Making the substitution that h = H/W and rearranging yields the follow- 
ing relation 


Hf H i'^ 


W„ 


w. 

1 


(37)' 


Since the heat exchange parameter q is defined as 


q 



3 


(38) 
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equation (37) can be rearranged to yield 


H f H i 

f 1 


(39) 


Since the values of W^, IL, and q are known, equation (39) can be 
solved by trial and error to obtain T^. The method is to assume a 
value for T^, find and H f by the procedure outlined in part 1, 

■ < . t *• . ■ *. f 

and compare the two sides of equation (39) • This method is repeated 1 
until covergence or satisfactory agreement is reached. 


3. Equations for Isentropic Flow From a Specified Chamber 
Condition to a Specified Pressure in the Nozzle 

With the temperature and pressure of the chamber known, the 
entropy, enthalpy, and other thermodynamic properties of the dissociated 
hydrogen gas in the chamber can be calculated by the procedure outlined 
in part 1. For isentropic flow between the chamber after heat removal 
and a specified point in the nozzle, y,. the change in the entropy per 
unit mass is zero; that is, 



or 




(4o) 


(4l) 
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Assuming that the pressure is known at position y, the temperature T , 

y 

can be found by a trial and error method. This involves making an 
estimate for T , performing the calculations in part 1, and then substi 

y 

tuting the results into equation (4l). This process is repeated until 
the entropies converge or until satisfactory agreement is reached. 


Once the temperature at y has been determined, the composition and 
other thermodynamic properties of the gas mixture are found by using 
the method outlined in part 1 involving equations (26) through (34). 
The velocity V may now be calculated using 

y 


Ah = -A 



(42) 


where use has been made of assumptions (l) through (5) listed in chap- 
ter I. Since h = H/W and = 0, equation (42) may be written 



V 

JL- 

2 


(43) 


or 




(44) 


Since the thermodynamic properties of the gas at f and y have been 
determined as outlined above, equation (44) may be used to calculate V 

y 
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Thus, by using the procedure outlined in this part, it is 
possible to find all of the desired conditions at point y providing 
the temperature and pressure in the reservoir (i.e., and P^) and the 
pressure at point y are known. 

Instead of designating the area ratio as the parameter 
defining the point of interest in the nozzle (chapter II), the pressure 
P was chosen as the designating factor. The calculation of the area 

y 

ratio at position y in terms of the pressure P will be shown in 

y 

part 5* 


4. Equations for Heat Addition in the Nozzle 
The problem in this segment of the nozzle is to find the 
thermodynamic properties of the gas mixture behind the heat addition 
plane. The conditions existing immediately in front of the heat addi- 
tion plane are known from the procedure described earlier in this 
chapter (part 3), and the amount of thermal energy being added to the 
gas flow is the same as that subtracted from the chamber (part 2). 

Across the heat addition plane, the mass, momentum, and energy conserva- 
tion equations must .hold. Penner (ref. 2) combines these conservation 
equations and derives a relation between the initial pressure and final 
pressure and a relation between the initial and final enthalpies. These 
relations are 


)l/2 


( T 2 + 1 ) 


% - -) 

~~t7 


- 4y 


2 T_ 


( 45 ) 


2y 


2 T, 
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and 



(46) 


Equations (45) and (46), along with equation (34), are used in the 
iteration scheme to provide better values for y^ as the iteration 
progresses . 

In order to start the iteration, it is necessary to estimate a 
value for This first estimate can be obtained by solving the 

Rayleigh relations as was done in chapter II, equation 17. 


and 



i 1 * Y i M i) 2 



(Vr) 


T 2 _ 2 M FE ( Y FE + - 1 ) 

t t " i 1 * ^ 



Dividing equation (47) "by equation (48) results in 


(48) 


1 

C\J 


1 1 + y 4si 

2 1 

1 

T 

T i 

m i j 

| 1 + y M 2 

i 

\ 

i + 1 ^4e 


( 49 ) 
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assuming that y = y = y as was assumed in chapter II. Since 1VL , 

X riL _L 

T , and T are known, equation (49) can he 'solved for the first 
1 2 

estimate of the Mach number behind the heat addition plane (M^). Once 
M^ had been calculated, the first estimate of the temperature behind 
the heat addition plane (Tpg) can be found by using 


T 

PE 


1 + 




( 50 ) 


where the stagnation temperature is, in reality, the initial 

2 


reservoir temperature' T\. 

The procedure to be followed is to calculate the first esti- 
mated value of the temperature T_ by using equations (4-7) through 

r ijj 

( 50 ). Then assuming that the first estimate for y^ equals y^, 
equation (45) can be solved with T^ substituted in place of Tg for 
the ratio Pg/P^. Next, using the first estimates for Pg and Tg, 
an improved value of y^ can be found from equation (34). With the 
improved estimates for Pg and y^, equation (46) can be solved for 
hg and, thus, Tg. With new values for Tg and. Yg, equation (45) 
again yields an improved value for P . The iteration is continued 
until the temperature, pressure, and ratio of specific heats have 
converged to definite values. After convergence the thermodynamic 
properties of the gas mixture immediately behind the heat addition plane 
can be found by the procedure developed in part 1. 
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5. Calculations of the Area Ratio 

Given the pressure, temperature, gas velocity, and all other 
flow properties at a point in the nozzle, the only remaining desired 
parameter of interest is the area ratio. This ratio can he calculated 
hy using the continuity equation and table II of NASA SP-3002 (ref. 6). 
From this table for a specified temperature and pressure, the compressi- 
bility factor Z and the ratio of the equilibrium speed of sound to 
the equilibrium speed of sound at 0° C can be found. The equation of 
state yields the density 


_P_ 
p ZRT 

and the ratio a/a Q yields the equilibrium speed of sound 


(51) 


a 



(52) 


Assuming that the Mach number is equal to unity at the throat, the gas 
velocity is given by the equilibrium speed of sound. NASA TN-D 1^5^ 

(ref. 7) reaches the conclusion that an excellent estimate for the 

* 

throat pressure P for equilibrium flow is 

P c ( y c * 1 I Y c - 1 
p ' 


2 


5 


(53) 
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where y q is the isentropic exponent in the combustion chamber and 
P c is the pressure in the combustion chamber. The throat temperature 
T can be calculated from 


T 

c 

* 

T 


Y +1 

c 

2 


( 5*0 


where !P is the chamber pressure. 

Using the continuity equation, the area ratio may be written 


as 



( 55 ) 


The density at the throat p is calculated by using equations (5l), 
(53), and (5^)« The density p and the velocity V are found by 
using equation (5l) and the procedures described in part 3. V is the 
equilibrium speed of sound calculated by using equation (52) with the 
values of T and P determined from equations (53) and (5*0. These 
parameters, when substituted into equation (55) 5 thus determine the area 
ratio at the point of interest. 

C. Example Check Calculations 

The method used to check the numbers calculated by the computer 
program was to compare the variables calculated during each step to 
values of the variables for the same initial conditions contained in 



NASA TN D-275 (ref. 8). The conditions used for one set of comparisons 
were: 


ex 


= 6b. 99 


■£ = 27.77 


A 

T. = 1+000° K 

1 

P^ = 30 atm 
q. = 0.25 

The comparisons for the above conditions are presented in the following 


tables . 
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COMPARISON OF NASA TN D-275 AND CALCULATED VALUES 
FOR SELECTED VARIABLES. . 

[From ref. 8] 


Variable 

NASA TN D-275 

Calculated 

1 Difference , percent 

Initial reservoir conditions 

T 

Uooo.o 

4000.0 

N/A 

P 

30.0 

30.0 

N/A 

h 

23911.1 

23979.6 

0.29 

m 

1.763 

1.7632 

.00 

^2 

.74929 

.74920 _ 

.01 

s 

24.2382 

24.2399 

.01 

p* 

17.0890 

17.2058 

.68 

Ip* 

3751.0 

3750.9 

.00 

Reference exit conditions 

T 

1347.0 

1342 . 19 

0.36 

P 

.0300 

' .0296 

1.33 

h 

4725.1 

4708.67 

- .35 

m 

2.016 

2.016 

.00 

^2 

.99999 

.99998 

.00 

s 

24.2382 

24.2399 

.01 

M 

4.617 

4.61676 

.01 

I 

1292.1 

1294.9 

.22 

sp 



Immediately in front of heat addition plane 

T 

1821.0 

1814.66 

0.35 

P 

.1 

0.099 

1.00 

h 

6596.1 

6574.3 

.33 

m 

22.015' 

2.01456 

.02 

y 

1.3067 

1.3208 

1.08 

^2 

• 99866 

.998575 

.01 

S 

24.2382 

24.2399 

.01 

M 

3.841 

3-839 

.05 





















Mgw <fOCQ-<:p^‘xiHg 
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COMPARISON OF NASA TN D-275 AND CALCULATED VALUES 
FOR SELECTED VARIABLES - Concluded 
[From ref. 8] 


r— 

Variable 


NASA TN D-275 


Calculated 


Difference, percent 


Immediately behind heat addition plane 


2735.0 

2746.04 

0.40 

.22334 

. 22161 

.78 

14448.1 

14435.5 

.09 

1.875 

1.878 

.17 

1.145 

1.148 

.26 

26.7426 

26.7335 

.03 

5977.78 

5994.89 

.29 


Exit condition with heat addition 


2470.0 

2472.18 

0.09 

O.O767 

0 . 0761 

.78 

11559. 

11508. 

.44 

1.93667 

1-93577 

.05 

0.9213 

0.9204 

.10 

26.7425 

26.7335 

.03 

3.855 

3.7815 

1.91 

1420. 

i4o4. 

1.13 
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D. Results and Conclusions 

The results of applying the mathematical relations developed 
in section B and used in the formulation of the computer solution are 
presented in figures 19 through b2 for the two reacting gas cases inves- 
tigated. For the first case heat was added to the nozzle gas flow with 
no heat being withdrawn from the combustion chamber (figs. 19 through 
30) and, for the second case, heat was simultaneously withdrawn from the 
combustion chamber and transferred to the nozzle gas flow (figs. 31 
through h2) . The ratios of velocity, pressure, specific impulse, tem- 
perature, Mach number, and molecular hydrogen mole fraction were plotted 
against the area ratio. The area ratio of the heat addition plane was 
varied by moving the plane- from the throat to the nozzle exit (fig. l4). 

For the flow of a dissociating diatomic gas, with heat addition 
in the nozzle but no heat withdrawn from the combustion chamber, the 
exit velocity ratio decreased as the plane of heat addition was moved 
from near the throat to the nozzle exit (fig. 19)* The decrease was 
small, however, with the velocity ratio remaining above 1.0 and with no 
"crossover" point, such as that found by Smith (ref. l) for external 
heat addition to the nozzle gas flow for an inert gas. Increased exit 
area ratio resulted in an increase in the exit velocity ratio (fig. 20) 
and also in a smaller net decrease when the heat- addition plane was 
moved from the throat to the nozzle exit. The exit pressure ratio in- 
creased as the plane of heat addition was moved from near the throat to 
the nozzle exit (fig. 21). Increased nozzle exit area ratio caused a 
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slight increase in the exit pressure ratio (fig. 22). The specific im- 
pulse ratio behaves similarly to the exit velocity ratio, thus the de- 
scription of the exit velocity ratio is applicable to the specific impulse 
ratio (figs. 23 and 2k ) . The exit temperature ratio behaved similarly 
to the exit pressure ratio (figs. 25 and 26). The ratio of the exit Mach 
number is always below 1.0 and decreases as the plane of heat addition is 
moved from near the throat to the nozzle exit (fig. 27). Increasing the 
exit area ratio causes a larger decrease in the exit Mach number ratio 
(fig. 28). The exit molecular hydrogen mole fraction ratio is below 1.0 
and decreases as the plane of heat addition is moved from near the throat 
to the nozzle exit (fig. 29). Increasing the exit area ratio causes an 
increase in the exit molecular hydrogen mole fraction ratio (fig. 30). 

Thus when heat is added to the flow in the nozzle, but not removed from 
the chamber, the dissociating diatomic gas shows an increase in ; exit 
velocity, exit pressure, exit temperature, and specific impulse. 

Comparison is made of the results of external heat addition 
to the nozzle gas flow for a dissociating diatomic gas (fig. 19 through 
30) with the results of M. L. Smith's analysis (ref. l) for external heat 
addition to the nozzle gas flow for an inert gas. The comparison shows 
similarities as well as differences. The tendency of the respective 
curves is the same. For example, if the analysis for the inert gas shows 
that a particular ratio decreases as the plane of heat addition is moved 
from near the throat to the nozzle exit, the same ratio will also decrease 
for the dissociating diatomic gas. The main difference is that the exit 
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velocity ratio for the dissociating diatomic gas case does not have the 
"crossover" that was shown by M. L. Smith in his analysis for an inert 
gas. The primary reason that the exit velocity ratio for the reacting 

gas remains above 1.0 is the fact that the chemical reaction of recombi- 

nation, causes larger amounts of heat to be added to the flowing gas up- 
stream of Smith* s "critical" point than are added downstream, thus 
causing a higher exit velocity. 

Figures 31 through k2 present the results for the case of a 
dissociating diatomic gas in which heat was transferred from the com- 
bustion chamber, by some means unspecified in this analysis, to the flow 
in the nozzle. The exit velocity ratio, which is always less than 1.0, 
was found to decrease as the plane of heat addition was moved from near 

the throat to the nozzle exit (fig. 3l). Increased exit area ratio 

caused an increase in the exit velocity ratio (fig. 32). The exit pres- 
sure ratio remained above 1.0 and increased as the plane of heat addition 
was moved from near the throat to the nozzle exit (fig. 33). Increased 
exit area ratio caused an increase in the pressure ratio (fig. 3*0 • The 
specific impulse ratio behaved similarly to the exit velocity ratio 
(fig. 35)* Increased area ratio caused an increase in specific impulse 
ratio (fig. 36). The exit temperature ratio remained above 1.0 ^ 
and increased as the plane of heat addition was moved from near the 
throat to the nozzle exit (fig. 37)* Increased exit area ratio caused 
an increase in the temperature ratio (fig. 38). The exit Mach number 

ratio decreased as the plane of heat addition was moved from near the 
throat to the nozzle exit and remained below 1.0 (fig. 39)* Increased 
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exit area ratio caused a decrease in the exit Mach number ratio (fig. 4o). 
The exit molecular hydrogen mole fraction ratio behaved similarly to the 
exit velocity ratio (figs, hi and 42). 

Comparison of the, results of the analysis of heat transfer from 
the combustion chamber with the nozzle gas flow for both the inert gas 
case and dissociating diatomic gas shows several similarities as well as 
a few differences. Comparison of the two cases shows that the tendency 
of the curves is the same. For example, the exit velocity ratio in both 
cases decreased from. a maximum close to 1.0 at the throat to the lowest 
value at the nozzle exit. The same similarity holds true for the re- 
spective curves of both cases. A direct comparison of the. inert gas. and 
dissociating diatomic gas (as far as the specific impulse ratio) shows 
that the specific impulse ratio for an inert gas varies from 0.94 to 
0.84 while the reacting gas remains close to 1.0 (fig. 43). Results 
of the inert gas analysis show that transferring heat from the com- 
bustion chamber to the nozzle causes a decrease in, both exit velocity 
and specific impulse but results in an increase in exit pressure and 
exit temperature. The, dissociating diatomic gas behaves in a manner 
similar to the inert gas, but the magnitude of the decrease is not as 
great. The damping effect is caused by the recombination reactions and 
associated heat release in the diatomic gas. 
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Figure 3.- Flow diagram for inert gas 


case 



Figure 3-- Continued 
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Figure 3*- Continued. 
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MA==0 

IX- 

MA=MA+1 


DELT=T-S(MA,l) 



MB=MA-1 

PAM=DELT/(S(MA,l)-S(MB,l) ) 

AKP= S ( MA , 2 ) + ( S ( MA , 2 ) - S ( MB , 2 ) ) * PAM 
SHR=S(MA,3)+(S(MA,3)-S(MB,3) )*PAM 
SH2=S(MA,4)+(S(MA,10-S(MB,4) )*PAM 
HHR=S (MA, 5 )+(S(MA, 5 )-S (MB, 5 ) )*PAM 
HH2=S(MA,6)+(S(MA,6)-S(MB,6) )*PAM 
CPH=S(MA,T)+(S(MA,T)-S(MB,7))*PAM 
CH2=S(MA,8)+(S(MA,8)-S(MB,8) )*PAM 

BETA= ( AKP/ ( U*P+AKP ) ) **0 . 5 

XH=2 . *BETA/ ( 1 . +BETA ) 

XH2= ( 1 . -BETA ) / ( 1 . +BETA ) 

W=WH*XH+WH2T*XH2 
C =CPH*XH+CH2*XH2 

T=C p /(C -H) 

H=XH* ( HHR+DHFO ) +XH2*HH2 
H=H/W 

S=XH* ( SHR-R* £n ( XH ) )+XH2*(SH2-R*5,n(XH2) )-R*AnP 
S=S/W 



Figure l6.- Flow diagram for Subroutine Scan. 











KEY=1 

INAN=0 



Figure 18.- Flow diagram for 


Imi 


END 

STOP 


Dissoc 























TC OMP=TCOMP+ZA ( MOT ) 



CALL SCAK 

IK: 

TCOMP, PR 

OUT: 

XH2C0 , WBCO, GAMCO 


HBACO, EKTCO 


HBACO : HBCOMP 


TCOMP=TCOMP-ZA(MOT) 

MOT=MOT+l 


TSTAR=2 . 0*TC0MP/( 1+GAMCO ) 

/GAMCO \ 

PSTAR=PR/ ( GAMpi J VGAMCO-lj 



CALL COMPUT 

IN: 

TSTAR, PSTAR, KEY 

OUT: 

RSTAR, ASTAR 


Figure l8.— Continued 
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TADD=TADD+ZA ( KYZP ) 

~r~ 

CALL SCAN 


IN: TADD, PADD 

OUT: XH2ADD, ¥ADD, GAMADD 

HBADD, ENTADD 



CALL COMPUT 

j IN: TADD, PADD, KEY 

OUT: RADD, AADDT 

l 

VADD = GRAV*(2.0* ( HBCOMP-HBADD ) ) 
ARADD = RSTAR*ASTAR/(VADD*RADD) 



Figure l8 . - Continued . 
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Figure 18.- Continued. 
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Figure l8.- Continued. 






TINT=TKONT 


GPAEE=GAMINT 
PPARE=PINT 
TPARE=TINT 
m3M=NUM + 1 


© 


PNEWR=0.0 , 
MOB= ' 3 




Figure 18.- Continued. 




CALL SCAN 


IN : TR , PNEWR 

OUT: XH2NEW, WHNEW, GAMNEW 

HBANEW, ENTNEW 



Figure 18.- Continued. 
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CALL SCAN 


IN: TFIND, PFIND 

OUT: XH2IND, WFIND, GAFIND 

HBFIND, ENFIND- 



Figure 18.- Continued. 
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Figure 18.- Continued. 
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Figure 18.- Concluded. 
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Figure 23-- Specific 
values of the heat 
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impulse ratio versus area ratio of heat addition for various 
exchange parameter (reacting gas with heat addition). 
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Figure 2 . 6 .- Exit temperature ratio versus area ratio of heat addition for three values 
of the exit area ratio (reacting gas with heat addition). 
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